An exact relation between Eulerian and Lagrangian velocity increment statistics 
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We present a formal connection between Lagrangian and Eulerian velocity increment distributions 
which is applicable to a wide range of turbulent systems ranging from turbulence in incompressible 
fluids to magnetohydrodynamic turbulence. For the case of the inverse cascade regime of two- 
dimensional turbulence we numerically estimate the transition probabilities involved in this connec- 
tion. In this context we are able to directly identify the processes leading to strongly non-Gaussian 
statistics for the Lagrangian velocity increments. 
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Introduction The relation between Eulerian and La- 
grangian statistical quantities is a fundamental question 
in turbulence research. It is of crucial interest for the un- 
derstanding and modeling of transport and mixing pro- 
cesses in a broad range of research fields spanning from 
cloud formation in atmospheric physics over the disper- 
sion of microorganisms in oceans to research on combus- 
tion processes and the understanding of heat transport 
in fusion plasmas. The recent possibility to assess the 
statistics of Lagrangian velocity increments by experi- 
mental means has stimulated investigations of rela- 
tions between the Eulerian and the Lagrangian two-point 
velocity statistics. Especially, the emergence of intermit- 
tency (i.e. the anomalous scalingof the moments of the 
velocity increment distributions [3[) in both descriptions 
and its interrelationship is of great importance for our 
understanding of the spatio-temporal patterns underly- 
ing turbulence. A first attempt to relate Eulerian and 
Lagrangian statistics has been undertaken by Corrsin [3| , 
who investigated Eulerian and Lagrangian velocity cor- 
relation functions. Recently, the Corrsin approximation 
has been reconsidered by [5|, where it has become ev- 
ident that it is a too crude approximation and cannot 
deal with the question of the connection between Eule- 
rian and Lagrangian intermittency. Further approaches 
to characterize Lagrangian velocity increment statistics 
are based on multifractal models [a, 0] which also have 
been extended to the dissipation range In Q a di- 
rect translation of the Eulerian multifractal statistics to 
the Lagrangian picture is presented. In this approach 
a non-intermittent Eulerian velocity field cannot lead to 
Lagrangian intermittency. This statement is in contrast 
with the experimental results of Rivera [§], as well as 
numerical calculations performed for 2d turbulence [lo| . 
Motivated by this fact we have derived an exact relation 
between the Eulerian and the Lagrangian velocity incre- 
ment distributions, which allows to study the emergence 
of Lagrangian intermittency from a statistical point of 
view. 

Connecting the increment PDFs The quantities of inter- 



est are the Eulerian velocity increments 
Ue = v{y + x,t) ~v{y,t), 



(1) 



where the velocity difference is measured at the time t 
between two points that are separated by the distance x, 
and the Lagrangian velocity increment 



ui 



v(y + x{y, T, t), t) - v{y, t ~ r). 



(2) 



In the latter case the velocity difference is measured be- 
tween two points connected by the distance x{y,T,t) 
traveled by a tracer particle during the time interval r. 
In both cases v is defined as the projection t; • of the 
velocity vector on one of the axes (i = x,y, z in 3d and 



« = a:, y in 2d ) of the coordinate system (see e.g. |12l.ll3l|). 
In the case of an isotropic fiow the results do not depend 
on the chosen axis, however we do not have to make this 
assumption yet. Additionally, we define the velocity in- 
crement 



Uei = v{y + x{y, T, t),t) - v{y, t), 



(3) 



which is a mixed Eulerian-Lagrangian quantity because 
the points are separated by x but the velocities are mea- 
sured at the same time. The properties of this quantity 
have been investigated in [ij]- Finally, we introduce 



v{y,t) - v{y,t- r). 



(4) 



measuring the velocity difference over the time r at the 
starting point of the tracer. Following 15, 1^ we define 
the so called fine-grained PDF for Ue as 



fe{Ve;X,y,t) = S{Ue - We). 



(5) 



where Ue is the random variable and is the independent 
sample-space variable. The fine-grained PDF describes 
the elementary event of finding the value Ve given the 
measured Ue — v{y -\- x,t) — v{y,t). The relation to the 
PDF is determined by 



fe{Ve]X,y,t) = {fe{Ve\X,y,t)), 



(6) 



2 



where the brackets denote ensemble averaging. The 
quantity fe is a function with respect to the variables 
x,y,t and a PDF with respect to the variable v^. In 
analogy to ([5]) we can define fine-grained PDFs for all 
other quantities. Now we want to derive an exact rela- 
tion between the fine-grained PDFs fe and fi. This task 
can be split up into two steps. First we have to replace 
the distance x by the trajectory x{y, t, t) of a tracer in 
order to translate from Ue to Uei- This is done by 

fel{Ve;y,T, t) 

= J dx S{x{y,T,t) ~ x)feive;x,y,t). (7) 

We see that during this operation the sample-space vari- 
able is not affected. The subscript in/ej denotes the fact 
that we now have Ve = Uei instead of Ve — u^. In 
the second step we have to connect /; and fei- From 
the definitions of the increments ([2]) - ([1]) we see that 
ui = Uei + Up and therefore the fine-grained PDF for ui 
is given by the fine-grained distribution of the sum of u^i 
and Up. To that end we have to multiply fei{ve]y,T,t) 
with fi,i{vi,;y,T,t) = 5{up — Vp) to get the fine-grained 
joint probability of finding u^i and Up at the same time. 
Subsequent application of J dve J dvp d{vi — {ve + Vp)) 
leads to 

fi{vi]y,T,t) 

= j dVe fpivi ~ Ve;y,T,t)fel{Ve;y,T,t). (8) 

To derive the corresponding PDFs we have to perform 
the ensemble average. In case of equation ([5]) we obtain 

fi{vi;y,T, t) 

^{JduJpiu,-u,,y,rM.iv,,y,r,t)) 

= j dVe fp{vi - Ve\Ve;y,T,t)fel{Ve;y,T,t). (9) 

In the last line we used the general relation p{a, b) — 
p{a\b)p(b) valid for two random variables in order to ex- 
tract fei from the average. We can treat ([71) in a similar 
manner. This leads us to 

fel{Ve;y,T,t) 

= (^J dx 6{x{y,T,t) - x)fe{ve;x,y,t)^ 

= j dx {5{x{y,T,t) ~ x)\ve)fe{ve;x,y,t). (10) 

Inserting equation ^ into pO|) shows that the Eulerian 
and the Lagrangian velocity increment PDFs are con- 
nected via the transition probabilities pa — {S{x{y, r, t) — 
x)\ve) and = fp{vi - Ve\ve;y,T,t) = fp{Vp\Ve;y,T,t). 




FIG. 1: The plot shows Pa (r I lie; T") (top) and Pa{vp\ve; t) (bot- 
tom) for r = 0.09Ti. 

Before we connect both equations we introduce some sim- 
plifications. In most experiments and numerical simula- 
tions dealing with Lagrangian statistics the flow is as- 
sumed to be stationary and homogeneous. In this case 
we may average with respect to y and t and hence the 
dependence on this parameters in (jlOp and ([9]) drops. 
Under the assumption of isotropy fei depends only on 
r = \x\. Therefore, we can introduce spherical coordi- 
nates in (fTO|l and integrate with respect to the angles. 
As mentioned before, in the case of isotropy the statisti- 
cal quantities do not depend on the chosen axis. Finally 
we arrive at 

= J dVePb{vi - Ve\Ve;T) J dr p a{r\v e] t) f e{v e] r) . 

/e,(i'.:r) 

(11) 

For convenience, we included the integrated functional 
determinant (27rr in two and Airr"^ in three dimensions ) 
in Pa- In this case Pa{r\ve', t) is a measure for the turbu- 
lent transport and gives the probability of flnding a tracer 
traveling the absolute distance r within the time interval 
T. Multiplication by fe{ve',r) and subsequent integration 
over the whole r-range mixes the Eulerian statistics from 
different length scales r weighted by Pa to form the PDF 
fei{ve',T) for a fixed time-delay t. The occurence of the 
condition in pa shows that this weighting depends on Ue- 
changes by Up. The transition probability pi, incorporates 
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FIG. 2: The figure shows the impact of the transition probabilities on the Eulerian PDF. The left part of the figure shows 
fe{ve',r) for r — 0.06,0.12,0.3. These PDFs are transformed into f{ve;r) (middle) by the transition PDF presented in the 
upper half of Fig. [T] Subsequently f{ve;r) is converted into fi{vi;T) (right picture) by the second transition PDF from Fig. 
[l] In both cases r = 0.09T/. In plot II and III the points denote the reconstructed PDFs based on (|10() and ^ and the lines 
denote the PDFs directly computed from the Lagrangian data. 



the fact that during the motion of the tracer particle 
the velocity at the starting point Multiplying fe.i(ve',T) 
with pb and integrating over Ve sorts all events where 
Uei +Up — ui into the corresponding bin of the histogram 
for ui . We want to stress that equations (fTO|) and ([9]) and 
except from symmetry considerations also equation pT|) 
are of a purely statistical nature and, as a consequence, 
hold for quite different turbulent fields. They are valid 
for two-dimensional as well as three-dimensional incom- 
pressible turbulence but can also be applied to magnc- 
tohydrodynamic turbulence. The differences in the de- 
tails of these turbulent systems, which are connected to 
the presence of different types of coherent structures like 
localized vortices in the case of incompressible fluid tur- 
bulence or sheet-like structures in magnetohydrodynamic 
turbulence, are therefore closely related to the functional 
form of Pa and pb ■ 

Two-dimensional turbulence In this section we want to 
estimate numerically the two transition PDFs in (fTT|) for 
the case of the inverse energy cascade of two-dimensional 
turbulence. The data are taken from a pseudospectral 
simulation of the inverse energy cascade in a periodic 
box with box- length 27r Recapitulating the deriva- 
tion of (fTT|) we see that we need the velocity at the start 
and the end point of a tracer trajectory at the same time 
to estimate the transition PDFs. Therefore we have to 
record the velocity at the starting points of the tracers 
additionally to their current position and their current 
velocity. In Fig. [T] the transition probabilities Pa and pb 
are depicted for r = 0.09Tj, where Tj denotes the La- 
grangian integral time scale Uj. We have chosen this 



rather small time lag as an example because in this case 
the deviation of the Lagrangian increment PDF from a 
Gaussian is significant. The transition probability can 
be approximated by 

Pa{r\ve;T) = N{ve,T)rexp[~{r - m(we, r))^/CT^(t;e, r)]. 

(12) 

For small Ve we have m{ve, t) ^ a(r)|ue | and a'^{ve, t)) ^ 



I3{t). From the functional form of pa one can see that 
the transport of the tracer particles is of probabilistic 
nature. For any fixed Ve the tracers travel different dis- 
tances during the same time r. We also see a strong 
dependence on the condition which can be interpreted 
as deterministic part of the turbulent transport. This 
distinguishes it from pure diffusion and directly shows 
that the widely used Corrsin approximation is violated. 
In the case of deterministic transport pa would be pro- 
portinal to 5{r — m{ve,T)). A good approximation for pb 
is given by 

Pb{vp\ve;T) = 7V(i;e,r)exp[-(up - m(t;e, r))^/cr^(we, r)] 

(13) 

with m{ve,T) — q;(t) tanh(/3(T)t;e) and a'^{ve,T) = 
7(t)(1 + 5{T)\ve\)- For small Ve we see a strong nega- 
tive correlation between Vf, and Vp (here i&nh.{ve) ~ Ve)- 
Both quantities tend to cancel in this case. This neg- 
ative correlation between the sample-space variables in 
Pb is connected with the sweeping effect. For a tracer 
starting with vi travelling the time r without changing 
its velocity we have u^i = f i — f 2 when the velocity at 
the starting point changes during r from vi to V2- In 
this case we have Up = V2 — vi = —u^i. This corre- 
sponds to an idealized situation but it gives a hint at the 
cause of the negative correlations in pb- For larger Ve the 
correlation decreases. This is captured by the fact that 
tanh(/3(r)?;e) const for large Ve- For both transition 
PDFs we observe that for increasing r the dependence 
on their conditions vanishes. 

Now we want to turn to the question how the transi- 
tion PDFs transform the Eulerian PDF fe{ve',r) into the 
Lagrangian PDF fi{ui;T). This process is depicted in 
Fig. [2l The left part of the figure shows several examples 
of fe(ve]r). Applying J dr pa{r\ve',T) (see equation [TT|) 
superposes different Eulerian PDFs fe{ve',r) with differ- 
ent variances leading to the triangular shape in the semi- 
logarithmic plot of the new PDF fei{Vf.\T) (middle of 
Figl2]). During the transition from fei{Vf.\T) to fi{vi;T) 
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FIG. 3: The figure shows Pa{r;T) for t = 3.5r^, 14t^, 28r^ 
for three-dimensional turbulence. In the inset Pa is depicted 
for the two-dimensional case with r = 0.22Tr, 0AATi,0.9Ti 



translation rule in 0] we would have a fixed relationship 
1^1 VeT (Sur correspouds to in our notation) between 
the time lag r and the distance traveled by a tracer parti- 
cle during this time. This would correspond to choosini 
Pa S(r — Vst) in pT|) . The additional assumption 
Svt 6ur {ui ^ Ue in our notation) that the velocity 
fluctuations on the time scale r are proportional to the 
fluctuations on the length scale r could be incorporated 
in our framework by choosing pf, ^ 5{vi — Ve) leading to 

fl[vi]T) = J dveS{vi - Ve) J dr 5 [r - V ^t) f e{v e] r) 
= fe{vi;viT). (14) 

The exponents for the Lagrangian structure function ex- 
ponents can now be obtained by making use of the Mellin 
transform 



the variables Vp and Ve are added to form vi . The previ- 
ously described observation that Vp and Ve tend to cancel 
each other for small Ve leads to a stronger weighting of 
very small vi so that the new PDF /; [vi ; t) is strongly 
peaked around zero (right part of Figl2l). In contrast to 
the center of the distribution the tails seem not to be 
influenced significantly by ph{vp\ve', t) . This is in agree- 
ment with the fact that for large Ve the correlation be- 
tween Ve and vi decreases. 

Three-dimensional turbulence To get an impression of the 
transition probabilities in three dimensional turbulence 
we used the data provided by 0, [13 to calculate the PDF 
Pa{r]T) — J dVf, pa{r\ve]T) for different r. The result is 
depicted in Fig. [3l We see that as in the two-dimensional 
case the Lagrangian time scale r is related to the Eu- 
lerian length scales by a PDF. This well known result 
shows that in principle it is not possible to connect them 
by a Kolmogorov type relation like r ^ r / 5ur [7]. In this 
relation r is a typical eddy turnover time connected to 
eddies of length scale r. In our example we have chosen 
three different values of the time delay r taken from the 
inertial range. Even when the maximum of pa is at a dis- 
tance r which lies in the Eulcrian inertial range there are 
significant contributions from very small and very large 
r. From this observation we can conclude that due to 
the turbulent transport the Lagrangian statistics is in- 
fluenced by contributions from the Eulcrian integral and 
dissipative length scales. 

Relation to the MuUifractal approach The characteriza- 
tion of the relation between Eulerian and Lagrangian 
PDFs with the help of the conditional probabilities 
Paiflve'TT) (Eulerian to semi-Lagrangian transition) and 
Pbivp\ve',T) (semi-Lagrangian to Lagrangian transition) 
allows one to recover a well known multifractal approach 
for relating Eulerian and Lagrangian structure function 
exponents 01 ■ As a side product, a simpler formula 
for the Lagrangian structure function exponents of this 
multifractal approach will be obtained. Regarding the 



feiVe,r) 



1 



dn Se{n)Ve " 



(15) 



with Se{n) = Ae(ri)r''=("^. Here, we use the same nota- 
tion as in Please note, that it will be not necces- 
sary to know the amplitudes Ae(n) but only the Eulerian 
scaling exponents Ce("-)- Using Eqn. (fT4|) and the Mellin 
transform we obtain 



fi{vi;T) = — 

'"I J -too 



dn Ae{n)T' 



C.(n)^Ce(«)-n^ 



(16) 



This Lagrangian PDF is now inserted into the inverse 
Mellin transform to obtain the Lagrangian structure 
functions 



Si{n) = 



dvi vlfi{vi;T) 



(17) 



dvi —Vi 

VI 



djAe{])T^^^'^v'l'^''^-'. 



Now we substitute = j - Ce(j), dj' = (1 - djCe{j))dj 
and denote the inverse function by j = ]{]')■ Thus we 
have 



Si{n) = 



' 1 
dvi —v^ 
ovi 



df Si{f)iSvi)-^' (18) 



with Si{j') 



AdlL 



TUQj^T^''^^'' and j' = j - Ce(i) and we 
obtain for the exponents 



0(" - CE{n)) = Ce{n) 



(19) 



It remains to show that this relation (|T9l) is identical to 
the formulas derived in Biferale et al. ^TJ. To see this, 
we shortly repeat the multrifractal approach which starts 
with the Eulerian structure function exponents 



Ce(p) =inf (p/i- 

h 



S - De{h)) ^ phi + S - Deihl) (20) 



5 



and p = D'^{h*). Here Df,{h) is the Eulerian singular- 
ity spectrum and /i* is the value where the infimum is 
achieved. The assumption r = v^t appears now in the 
denominator of the expression of the Lagrangian struc- 
ture function exponents 



^'^P^ = 'f[ rrn 



(21) 



From this it follows 

Clip-Ceip)) 

' {p - pK - 



inf 

h 



i + D,{hl))h + 2,-De{h) 
1 - h 



(22) 



In order to find the infimum we differentiate with respect 
to h 

D'^{h:)~D'^{h)+D,{h:)-D,{h)~D'^{h:)K+D'Mh = 

(23) 



From this it follows that h\ = h*^ and 



CUp - Ceip)) ^pK + 3- D.iK) = Ce{p) 



(24) 



which recovers (|19l) . A consequence of (ITO)) and is 
that for a self-similar Eulerian velocity field (as we can 
find it in the two dimensional inverse energy cascade) we 
should find self similar Lagrangian PDFs. As mentioned 
above this is in contradiction to recent experiments [o^ 
and our own numerical simulations [lo| . 

Conclusion and Outlook We presented a straight forward 
derivation of an exact relationship between Eulerian and 
Lagrangian velocity increment PDFs. For the example 
of two-dimensional forced turbulence we were able to ex- 
plain how it is possible to observe strongly non-Gaussian 
intermittent distributions for the Lagrangian velocity in- 
crements. The two mechanisms in this context are the 
turbulent transport of the tracers leading to the mixing 
of statistics from different length scales and the veloc- 
ity change at the starting point of the tracers leading 
to a further deformation of the increment PDF. In com- 
parison we analyzed data from simulations of three- di- 
mensional turbulence. Similar to the two-dimensional 
case we demonstrated that Lagrangian time- and Eu- 
lerian length-scales are connected via a transition PDF 
that varies with the time scale. We where also able to 
show that the well known multifractal model for the La- 
grangian structure functions ia a limiting case of the pre- 
sented translation rule. The next step is to estimate the 
transition probabilities for three-dimensional turbulence 



as well as magnetohydrodynamic turbulence in order to 
get a deeper understanding of the influence of the under- 
lying physical mechanisms, especially the presence of co- 
herent structures on the translation process. In this con- 
text the question why intermittency in the Lagrangian 
picture is stronger in magnetohydrodynamics than in 
fluid turbulence [l^, although the situation is reversed 
in the Eulerian picture, will be addressed. 
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